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Abstract 

I derive completely analytically the time evolution and final abundances of the light 
elements (up to 7 Be) formed in the big-bang nucleosynthesis. This highlights an interesting 
physics taking place during the formation of light elements in the early universe. 



1 Introduction 

The most occurrent chemical element in the universe is hydrogen. It constitutes nearly three 
quarter of all baryonic matter. The next mostly wide spread element Helium-^ , constitutes 
about 25%. The other light elements and the metals occur very rare. Very simple arguments 
lead to the conclusion that it is very unlikely that A He, deuterium (D) and other light elements 
could be burned in the stars (see, for instance, PQ,[2])- Therefore the only sensible explanation 
of their abundance is that they were produced in the very early universe. It is clear that the 
essential amount of the helium could not be formed before the temperature dropped below its 
binding energy ~ 28 MeV and one can expect that the Big-Bang nucleosynthesis (BBN) took 
place when the temperature was not very different from ~ MeV, that is, somewhere in between 
seconds and minutes after the Big-Bang. Therefore Big-Bang nucleosynthesis (BBN), being based 
on the well understood physics, offers the possibility of reliable probe of the early universe (see, 
for instance, [I] , [2] , [H| , [Ij and references cited therein). The amount of the produced elements 
depend on the basic cosmological parameters and is very sensitive to the baryon density. The 
measured abundances combined with the CMB temperature fluctuation measurements provide 
us an unique opportunity to verify the reliability of the standard model of the universe evolution 

MM 

The element abundances, are usually calculated using computer codes (for instance, pub- 
lically available [Hj Wagoner code [HJ) and the abundances are presented as the function of the 
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baryon density. To understand the dependence of the element abundances on the cosmologi- 
cal parameters the semi-analytical and analytical description of BBN proved to be very useful 
[ZD];!!!]- in this paper I develop simple quasi-equilibrium analytical approach which allows to 
derive the final abundances of all light elements up to Beryllium-7 without using any computer 
codes. The accuracy of the results is very good for 4 ife, good for D and reasonably good for the 
other elements. I obtain analytical (not fitting) formulae describing the dependence of the abun- 
dances on the cosmological parameters and trace the time evolution of the element abundances 
before their freeze-out. This highlights an interesting and reach physics taking place during nu- 
cleosynthesis and allows to understand the physical reasons for the dependence of the abundances 
on parameters without practicing with computer codes. 



2 Freeze-out 

The amount of the produced helium depends on the availability of the neutrons at the time when 
the helium is formed. In turn, the neutron concentration is determined by the weak interactions 
which ensure the chemical equilibrium between the neutrons and protons at very early time. 
The weak interactions become inefficient when the temperature drops below few MeV. Around 
this time the neutrons chemically decouple from the protons and after that the ratio of their 
concentrations "freeze out" 1 . The nuclear reactions take place after that. Therefore, first we 
need to calculate the "freeze out" concentration of the neutrons. 

The main processes responsible for the chemical equilibrium between protons and neu- 
trons in the early universe are the weak interaction reactions: 

n + v ^ p + e~, n + e + ^-p + V. (1) 

Here v always means the electron-neutrino. To calculate the rate of these reactions one can use 
the Fermi theory according to which the matrix element characterizing 4-fermion interaction 
is equal: 

\M\ 2 = 16 (1 + 3g 2 A ) G F ( Pn ■ p v ) ( Pp ■ V e) , (2) 

where Gf = Tca w /y/2Mw ~ 1.17 x lO~ 5 /GeV 2 is the Fermi coupling constant, g& — 1-26 is the 
correction to the axial vector "weak charge" of the nucleon 2 and (p^ • pj) are the scalar products 
of appropriate 4- momenta entering the vertex 3 . Considering the process a + 6 — > c + g? of type 
we get the following expression for the differential cross-section of this interaction: 

2 / 2 \ ^' 

d(Tab = 1 \M\ ( jpc-Pd) -m 2 c m 2 d \ ^ 
dQ (8tt) 2 {p a + p b ) 2 \ (p a -Pbf - m\m\ J 
lr rhe above statement is, of course, literally true only if one neglects the neutron decay. 

2 This correction accounts for the possibility that the gluons binding quarks inside the nucleon can split into 
quark-antiquark pairs, which could give nonvanishing contribution to the weak coupling. 

3 The Fermi constant can be determined with a very good accuracy measuring the life time of the muon, while 
gA can be found only if one considers the interaction involving the nucleons. 
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where the integration over the phase space of c, d— particles has been performed. This expression 
is manifestly Lorentz-invariant and can be used in any coordinate frame. Note, that the 4- 
momenta of the produced particles are related to the 4-momenta of the colliding particles via the 
conservation laws: p c + Pd = Pa + Pb- Let is now consider the particular reaction n + v — > p + e~ 
at the temperatures around few MeV and below. In such a case the nucleons are nonrelativistic; 
hence 

{p n + Pvf - m\\ {p n ■ p u ) = m p e u ; 

(4) 



(P P ■ Pe) 2 ~ m 2 v m 2 e ~ m p \Jl - (m e /e e y 



m p e e v e 



where e v is the energy of the incoming neutrino and e e ~ e v + Q is the energy of the outgoing 
electron. The energy Q = 1.293 MeV, is released when the neutron "is converted" into proton. 
The formula is directly applicable only in empty space. However, at the temperatures above 
0.5 MeV there still present many e^— pairs and the possible final states for the electron are 
partially occupied. Because of the Pauli exclusion principle it reduces the appropriate cross- 
section by the factor (1 — n te ) = (1 + exp (— e e /T))~ l . Taking this into account and substituting 
(gj) into ©, © one gets: 



1 + 3(7 



a nv ~ ^—^AG 2 F e 2 e v e (1 + exp {-e e /T)y L (5) 



7T 



where we have neglected the chemical potential of the electrons. Note that the concentration of 
the nucleons is negligible compared to the concentration of the light particles at this time and 
therefore the spectrum of the light particles is practically not influenced by the above reactions. 
The nv— interactions taking place within time interval At in a given comoving volume, containing 
N n neutrons, reduce their total number by amount 



(6) 



where n €v = (1 + exp (—e u /T u )) , v v = 1 is the speed of the neutrinos and Ag ev is the phase 
volume element: 



Introducing the relative concentration of the neutrons 



Y - ^ n — Un (7\ 

n T\T I AT I ' V / 

N n + N p n n + n p 



and substituting (jSJ) in (jHJ) we finally obtain the following expression for the rate of change of the 
neutron concentration due to nv— processes 



dX n \ _ . v 1 + 3^1 



—X nu X n — ^ A G 2 F Q 5 J (1; oo) X n , 



dt ).. 2tt 3 
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where 




7(l:-x)= / <U,<fU,-\f \l-^' ,/Qr 



q 2 



1/2 




-1 








1 + e^ 9 



(9) 




and we have introduced the integration variable q = (e v /Q) + 1 = e e /Q. These expression is given 
in PJ. If we neglect the last multiplier into the integrand 4 and, taking into account that q > 1 
and (m e /Q) 2 ~ 0.15, expand the square root keeping only first two terms, the obtained integral 
can be calculated exactly and the result is 

, (1;oo) ,^ + ^ + ?C13) 11 _ ;| ^ ) w (10) 

where y = T v /Q. It is quite remarkable that this approximate expression reproduces the exact 
result with very high accuracy at all temperatures. For instance, at y > 1 the accuracy is about 
2%, improving to 1% and much better for y < 1. 

It is not difficult to understand why it is the case. Actually at low temperatures (y -C 1) 
this should be so since 5 exp (-Q/T) <C 1. On the other hand in the limit of very high temperatures 
(y ^> 1) the integral (jHJ) can be very well approximated if one neglects (m e /Q) and (Q/T u ) —terms; 
the result is 

77T 4 

J(l;oo)~^ 5 at ?/>l (11) 

Comparing this with the first term in (fTT)|) . which obviously dominates in this limit, we see 
that they coincide within 3%— accuracy since (45£ (5) /2) : (77r 4 /30) = 1.027. One can check 
numerically that in the intermediate range the accuracy of the approximate expression (fTU|) is 
better than 2%; for instance, at y = 0.7 it is about 1%. 

Substituting (fTUj) together with the numerical values of Gf, Q, expressed first in the 
Planck's units, into (jHJ) and then returning back to the usual units ([A] = sec -1 ) we infer that 

X nu ~ 1. 63y 3 (y + 0.25) 2 sec" 1 , (12) 

In this last expression the further simplifications were made. However, the reader can check 
himself that at all temperatures T u > 0.2 MeV its accuracy is never worse than 2 — 3% . Taking 
into account the experimental uncertainties in this accuracy looks very satisfactory. 



Similar by, we find that the rate of the reaction n + e + — > p + v is equal to 

1 + 3g\ 2 5 ( m e \ 
Ane = ~^T~ G fQ J I -oo; - — J , (13) 



4 This means that one ignores the Pauli's exclusion principle. 

5 We remind that before e ± — annihilation T = T„ and after that T = 1.4T„. 
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where J is the integral defined in (JUJ) with the limits of integration from — oo to — (m e /Q). If 
T u = T and T > m e , then X ne ~ \ nu . 

The rates of the inverse reactions: pe~ — > nv and pv — > ne + are related to the rate of the 
direct reactions (at T v = T) as 

Ap e = exp (-Q/T) X nu , \ pu = exp (-Q/T) \ ne , (14) 



"Freeze-out". The inverse reactions lead to the increase of the neutron concentration with 
the rate X p ^ n X p ; hence we can write the following balance equation for X n : 



~dT 



-K^ p X n + \_> n X p = -\^ p (l + er£) (X n - X?) (15) 



where \ n -+ p = A ne + X nu is the total rate of the direct reactions (JTJ) and X^ q = (1 + exp (Q/T))^ 1 . 
In deriving (JT5j) I took into account that the proton concentration X p = 1 — X n and used the 
relations (|T%|) assuming that T„ = T. 

The exact solution of this linear differential equation, with the initial condition X n — > X^ q 
as t — > 0, can be written in the following form 

X n (t) = X e n q (t) - J dtX? exp (- ^ X n ^ p (l + e'T^j di) (16) 

where dot denotes the derivative with respect to time. At small t the second term in this equation, 
characterizing the deviations from the equilibrium, is negligible compared to the first one. Inte- 
grating by parts, one gets that in this limit the solution (fTB^) can be rewritten as an asymptotic 
series in terms of the derivatives of X^ q 

Xn = X? | 1 - J- (l + e^)^ |g + .... ) (17) 



Therefore, if 



x e n q /x 



eq 



t 1 X n ->p, that is the rate of the reactions is very high compared to 
the inverse cosmological time, X n = X^ q , in complete agreement with the thermodynamical result. 
Much later, when the temperature significantly drops the "equilibrium concentration term" X^ 1 
goes to zero and at the same time the integral on the right hand side of (JTfij) approaches the finite 
limit. As a result the neutron concentration, instead of vanishing, as it would be in the case 
of chemical equilibrium, freeze-out at some value X* = X n (t — > oo) . 6 The freeze-out effectively 
occurs when the second term in (|17|) is of the order of the first one, that is, when the deviations 
from the equilibrium become significant. Assuming that this happens before e ± — annihilation and 

6 Note that if A would be decreasing not so fast, such that the integral in the exponent would diverge as t — > oo, 
then the overall integral term would also vanish in this limit. 
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after temperature drops below Q ~ 1.29 MeV (these assumptions can be checked a posteriori) 
one can put X n ^ p — 2X nu and neglect exp (-Q/T) in the obtained expressions. In this case the 



condition 



X eq /X eq 

n I n 



X n ^ p , defining the freeze-out temperature T*, takes the form 
^ 2 (^ + 0.25) 2 ~0.18k 1/2 



where = T*/Q. In deriving (fTS|) I used the formula (|T2*j) for X nu and took into account the 
relation between the temperature and cosmological time: 



'"(^f (t 1 ) ^ 1/2 ft (19) 



where k = ^ (db + §<?/) and gb, gf are the total numbers of the internal degrees of freedom, 
respectively, of all relativistic bosons and fermions. 

In the case of three types of neutrino (k ~ 3.54) ~ 0.65 and the freeze-out temperature 
is T* ~ 0.84 MeV. The equilibrium neutron concentration at this moment is X^ q (T*) ~ 0.18. 
Of course, this number gives only very rough idea about expected freeze-out concentration. One 
should not forget that at this moment the deviations from equilibrium are already very big and, 
in fact, X n (T*) exceed the equilibrium concentration at least twice. The most important thing 
which could be learned from this simple estimate is that the freeze-out temperature depends 
on the number of light species present in the universe at this time. Since T* oc k 1//8 , the more 
light species are present, the bigger is the freeze-out temperature and one can expect that more 
neutrons will survive after chemical decoupling from the protons 7 . In turn, later on nearly all 
these neutrons build 4 He; hence one can expect that if, for instance, in addition to known types 
of neutrino, there exist the other light particles, then the abundance of the primordial helium 
should be higher than in the case of three neutrinos. This can be easily understood if we take 
into account that the rate of the expansion of the universe (H = l/2i) at the given temperature 
increases if we have extra light particles (see ([19))): hence the freeze-out should occur earlier, 
when the neutron concentration is higher. For instance, in the extreme case of very big number 
of unknown light particles T* ^> Q and the expected concentration of the survived neutrons 
should be close to 50%, that is, there is one neutron per every proton. Later on these neutrons 
would bind the protons converting nearly all baryonic matter into 4 He. Of course, this would be 
in obvious conflict with the observational abundances of the light elements. Therefore, as we will 
see later, the primordial nucleosynthesis allows us to put rather strong restrictions of the number 
of light species. 



Now I calculate the freeze-out concentration more accurately. Since X^ q — > as T — > 0, 
this concentration is given by the integral term in (|Td^ where we have to take the limit t — ► oo. 
Changing the integration variable from t to y = T/Q (see ([19))) and taking into account that the 

7 If we use for freeze-out the simple criteria t ~ 1/A then we get T* oc k 1 / 6 , the result which is usually quoted 
in the literature. 
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main contribution to the integral comes at T > m e , when A n _> p ~ 2A n „ and X nu is given by (|12j) . 
we obtain 

• y » = r 2^1+ La/,)) ^ {- 5mk ' i/2 r * {x + 025)2 (i + e ^>) <2o) 

For the case of three neutrinos (k ~ 3.54) we get X* ~ 0.158. It is in a very good agreement 
with the results of more elaborated numerical calculations. The presence of extra light neutrino 
increases k by 2 ■ Anf ~ 0.58 and respectively the freeze-out concentration becomes X* ~ 0.163, 
Hence, two extra fermionic degrees of freedom (one for neutrino and one for antineutrino) lead 
to the increase of the freeze-out concentration by 0.5%. 

Neutron decay. In the above consideration I have neglected the instability of the neutron 
via decay 

n — > p + e~ + V (21) 

It was justified since the lifetime of free neutron r n = 885.7 ± 0.8 sec is rather large compared 
to the typical cosmological time at the moment of freeze-out (t* ~ 0(1) sec). However, later 
on the two-body reactions and inverse three-body reaction (}2Tj) become unimportant and the 
only remaining reaction reducing the amount of the neutrons is the neutron decay. As a result 
the relative concentration of the neutrons at t 3> t* is 

X n (t)=X*exp{-t/T n ) (22) 

Note that at late times one can neglect the degeneracy of the leptons which would increase the 
lifetime of the free neutrons; hence we can use the measured in the laboratory lifetime of the 
neutron quoted above. As we will see later the nucleosynthesis, as a result of which nearly all 
neutrons are captured in the nuclei, where they become stable, happens around t ~ 200 sec. It 
is a rather substantial fraction of the neutron lifetime and therefore the neutron decay changes 
significantly the amount of the survived neutrons and is important for the final 4 He— abundance. 



3 Deuterium bottleneck 

Complex nuclei are formed as a result of nuclear interactions of the baryons. For instance, 4 i?e 
could, in principle, be directly formed in many-body collisions; 2p + 2n ^ 4 He. However, the 
number densities at the time when this reaction can take place are too low and its rate is negligible 
compared to the rate of expansion. Hence, the light complex nuclei can be produced only in 
sequence of two-body reactions. The first step on this way is the deuterium (D) production: 

p + n ^D + 7 (23) 

There is no problem with this step since the rate of this reaction is very high and the "typical 
collision time" is, for sure, much smaller than the cosmological time (at t < 10 3 sec). Hence one 
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can expect that the deuterium should be in the local chemical equilibrium with nucleons. Let 
us define the deuterium abundance by weight as X D = 2n B /n Bt where n B is the total number 
of all baryons (nucleons) including those ones entering the complex nuclei. In the state of local 
chemical equilibrium the relation between X B and the abundances of the free neutron and protons 
(appropriately, X n = n n /riB and X p = n p /n B ) can be easily found with the help of the equilibrium 
Saha's formula (see, for instance, [2]): 



X D = 5.67 x 10^VoT A 3 / e V exp X p X n (24) 

where Bd = m p + m n — vtld ~ 2.23 MeV is the binding energy of the deuterium and the 
temperature is expressed in MeV. We have introduced here the normalized baryon-to-photon 
ratio 

77x0 = V10- 10 = /10- 10 , (25) 

which is related to the baryon contribution to the critical energy density fl B as 

tt B h 2 75 ~ 6.53 x 10~ 3 r7 10 . (26) 

where the Hubble constant h 75 is normalized on 75 km/ sec ■ Mpc. The abundance of deuterium 
at the temperatures about its binding energy is still very small. For instance, for T ~ 0.5 MeV, 
we get Xd ~ 2 x 10~ 13 . The reason for that is a very high entropy (number of photons) per 
baryon. Even at T <C B B there are still enough highly energetic photons with e > Bp which 
destroy the deuterium. Actually the number of these photons per one nuclei of the deuterium is 
about 

n 7 (e > B D ) BlTe-**/* _ 1q10 1 f Bd\ 2 ^ Bd/t (2?) 



n D n B X D Vw x d V T 

This number drops below unity at T < 0.06 MeV. Hence one can expect that the deuterium 
can be formed in significant amount only when the temperature is low enough, otherwise it is 
destroyed by the energetic photons. This also delays the formation of the other light elements 
as, for instance, 4 He. 

The binding energy of the helium-4 (28.3 MeV) is much higher that the binding energy 
of the deuterium; hence if helium would be in chemical equilibrium with neutrons and protons 
then one would expect that nearly all free neutrons would be captured in A He already at the 
temperature ~ 0.3 MeV. However in reality the helium abundance is still negligible at this tem- 
perature. This is because the rates of the reactions converting deuterium in more heavy elements 
is proportional to the deuterium concentration and is much smaller than the expansion rate until 
the deuterium abundance will increase and constitute the substantial fraction of the baryonic 
matter. Before that only the protons, neutrons and deuterium are in chemical equilibrium with 
each other. More heavy elements are decoupled and present in completely negligible amounts 
in spite of their high binding energies. This is known as "deuterium bottleneck". The size of 
the "bottleneck" which is proportional to X D should become big enough to allow the neutrons 
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and protons "to go through " and replenish the helium abundance in accordance with its high 
"equilibrium demand". Let us find when this happens. Using formula (J24|) we can express the 
temperature as a function of Xe>: 

°- 061 

TMeV {Xd) " (1 + 2.7X10-* In (XsAno)) (28) 

This relation is valid only when the deuterium is in chemical equilibrium with neutrons and 
protons, which as we will see is true until the moment when Xp reaches the value 1CT 2 . According 
to the formula (|28|) the deuterium abundance should change from 10 -5 to 1 when the temperature 
drops only in 1.5 times, namely, from 0.09 MeV to 0.06 MeV (for 7710 = 1). Therefore, the 
deuterium abundance should increase very abruptly around this time and one can expect that 
the nuclear reactions should become fast enough to proceed with formation of the light elements. 
The main processes converting the deuterium in more heavy elements are (see also Fig.l): 

1) D + D ^ 3 He + n, 2)D + D^T + p (29) 

The cross-sections of these reactions are known from experiments and the results are usually 
presented in terms of the effective rates vs. temperature 8 . In the temperature interval 0.06 + 0.09 
MeV these rates change not very much and we have 

(av) DD1 = (1.3 + 2.2) x 10- 17 cm 3 /sec, (ov) DD2 = (1.2 + 2) x 10" 17 cm 3 / sec (30) 

Considering the comoving volume containing No deuterium nuclei we find that the decrease of 
their number during the time interval At due to the reactions (}2*9"|) is equal to 

AN D = - (crv) DD n D N D At (31) 

Rewriting this equation in terms of the concentration by weight = 2Nd/Nb we obtain 

AX D = -h DD X 2 D At (32) 

where 

Add = (Mddi + Mdd 2 ) "b ~ 1-3 x 10 5 K (T) Tl IevVw sec- 1 (33) 

and K (T) is the numerical coefficient which changes from ~ 1 to ~ 0.6 when the temperature 
drops from 0.09 MeV to 0.06 MeV. It is clear that the substantial amount of the available 
deuterium can be converted into helium-3 and tritium within the cosmological time t only if 
AXd — (1/2)AddAd^ ~ Xd', hence the "deuterium bottleneck opens wide" only when the 
deuterium concentration reaches the value 

X{? ~ l 'l X lQ 5 ~ 1.5 x 10-Vo 1 (1 - 7 x 10- 2 ln7? 10 ) (34) 

VlO+MeV l A Dj 



3 The appropriate rates are cited in jllj . More recent data can be found on internet. 
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Deriving this formula I used the relations (|19|) with k ~ 1.11, and (|28j) ; the obtained equation 
was solved by iterations assuming that 1CT 1 < r] w < 10. 

After deuterium abundance reaches the value given by (JMj) everything proceeds very fast. 
In fact, if ?7io = 1 then according to (J28j) the equilibrium concentration should increase from 
10~ 4 to 10 -2 when the temperature drops from 0.08 MeV to 0.07 MeV. This increase of Xjy 
means that the reaction rates converting the deuterium to more heavy elements, which are pro- 
portional to Xp, at T ~ 0.07 MeV become 10 4 times bigger than the rate of the expansion. It is 
clear that this system is far from the equilibrium and the deuterium supplied by pn— reactions "is 
converted" very fast to more heavy elements. This doesn't allow the deuterium concentration to 
increase to the values bigger than 10~ 2 . The details of the nonequilibrium processes are described 
by a complicated system of kinetic equations which can be solved only numerically. In Fig. 2 the 
results of numerical calculations for the time evolution of the element abundances in the universe 
with VL B h% b ~5x 10~ 2 are shown 0. 

Below I present the calculations which explain the time behavior of these abundances and 
derive the formulae for the final freeze-out abundances of light elements up to 7 Be. This includes 
4 He, deuterium (D) , helium-3 ( 3 He) , tritium (T) , Lithium-7 { 7 Li) and Beryllium Be) . The 
other light elements as, for instance, 8 Li, 8 B etc. are produced in very small amounts and will 
be ignored. 

The most important nuclear reactions involving the light elements are schematically de- 
picted in Fig.l, which I recommend to keep in front of the eyes reading the rest of the paper. 




figure 1 



In this Figure to every element corresponds its own "reservoir" . All these "reservoirs" 
are connected by "one- way-pipes". Every "pipe" corresponds to an appropriate nuclear reaction. 
I write only the initial elements involved in the reaction, since the outcome can be easy inferred 
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from the picture. The "thickness of the pipe" through which the element a "escape from the 
reservoir" as a result of the reaction ab — > cd is proportional to rate of this reaction 

X a /X a = —A b l \ ab X b (35) 

where \ a b = ( (JV ) a b n B anc ^ ^ is the mass number of the element b; for instance, A = 4 for 
4 He and A = 7 for 7 Li, 7 Be. Of course, the appropriate "pipe" is efficient only if X a /X a > t~ l . 
As we have already seen the D— and p, n— reservoirs" are in equilibrium with each other and 
decoupled from the rest at the temperatures above 0.08 MeV ("deuterium bottleneck"). However 
when the temperature drops below 0.08 MeV the " DD— pipes open" and become very efficient 
in converting an extra deuterium supply from "np— reservoir" into more heavy element. Finally 
nearly all free neutrons disappear entering more heavy elements where they become stable. After 
that the concentrations of the elements in the appropriate "reservoirs" freeze-out and the "final 
abundances" survive. This is a general picture and now I proceed with detailed calculations and 
consider the formation of every element separately. 



4 Helium-4 

As soon as deuterium concentration increases to X { £ given by (|3*3j) the formation of the other 
light elements begins. This happens at the temperature (see 



T$ eV ~ 0.08 (1 + 7 x 10" 2 In 7710) (36) 

at the moment of time 9 

*«~ 206 (^_-J (1 -0.14 In ^o) (37) 

Of course, the nucleosynthesis does not happen instantaneously. Moreover at the beginning the 
rate of deuterium production in reaction, pn —* Dj, is substantially higher that the total rate of 
the deuterium "annihilation" in reactions (|29j) . namely, 

-r—YT - 10 " ("IF") 2 (38) 

where I used the experimental value for the ratio X pn /^DD, which is about 10~ 3 at TueV — 
0.07 4- 0.08 and put X n = 1 - X p ~ 0.16. 



As it follows from ()38|) before the deuterium concentration reaches its maximal value 
X D ~ 10~ 2 the deuterium production dominates over deuterium destruction and the deuterium 
abundance continues to follow its chemical equilibrium track given by (|24j) . According to (|28|) 



9 Note that rW and depend on the exact value of only logarithmically and therefore not very sensitive 



to the exact value of X D . 
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the concentration ~ 10 2 is reached very fast after namely, when the temperature drops 
from 0.08 MeV to 0.07 MeV (for r/ 10 = 1) , that is, with 



AT 

At ~ 2t (i) ^- ~ 50 sec 



(39) 



time delay after When this concentration is reached the two-body DD— deuterium destruction 
become more efficient than the pn— deuterium production and Xp begins to decrease 10 (see Fig. 2). 
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figure 2 
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The concentration of the free neutrons during this period strongly decreases and they go 
first to the "deuterium reservoir" and then proceed further "through the pipes" forming heavy 
elements. For most neutrons the "final destination" is the " 4 He— reservoir". 

Why it is so can be understood even without analyzing the rates of the intermediate 
reactions. Actually, if 4 He would be in the equilibrium with the other light elements it would 
be dominating at low temperatures because of its high binding energy (28.3 MeV) , which is 
four times bigger than the binding energies of the intermediate elements, 3 He (7.72 MeV) and 
T (6.92 MeV). The system which is away from equilibrium always tends there in a quickest 
possible way. Therefore, most of the free neutrons will be capture into A He— nuclei because its 
equilibrium demand is the highest. 

The reactions proceed in the following way. First, the deuterium is converted into 3 He 
and T in reactions (|29|). After that tritium interacts with deuterium and produce the helium-4 

10 The deuterium photo-destruction can be completely neglected after that. It is clear if we note that if there 
would be only photo-destruction processes alone then the deuterium concentration would continue to increase. 
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nuclei: 

T + D — > 4 He + n (40) 

As a result two neutrons out of three are captured into the 4 He— nuclei and one comes back into 
"rap— reservoir" . 

The 3 He— nuclei can interact either with free neutrons and then proceed to "T— reservoir" , 

3 He + n ^ T + p, (41) 
or with deuterium going directly to n4 He —reservoir 

3 He + D^ A He + p (42) 
The ratio of rates for these reactions is 

~ b ; (43) 

^HeD^HeXo Xjj 

hence at the beginning " 3 HeD— pipe" is inefficient compared to n3 Hen— pipe" and most of 3 He— 
nuclei are converted into tritium. Only when the concentration of the free neutrons drops below 
the deuterium concentration (which is always smaller that 10~ 2 ), the rate of the reaction (j42j) 
converting 3 He directly into 4 He becomes bigger that the rate of the reaction (j4*Tj) . It follows 
from here that most of the neutrons will go into A He— nuclei either along np — > D — > T — > 4 He 
or np — > D — > 3 He — > T He way. Finally, in about 50 4- 100 sec after the beginning of 
nucleosynthesis nearly all neutrons (with the exception of very small fraction < 10 -3 ), end up 
in 4 He— nuclei. Therefore, the final 4 He— abundance is completely determined by the amount of 
the available free neutrons at the time when DD-reactions become efficient, that is at t ~ 
Because half of the total weight of 4 He is due to the protons, its final abundance by weight should 
be 

X( He = 2X n (>«) = 2X1 exp {~) ( 44 ) 
Substituting here X* from (|2U|) and from (jH7|) we obtain: 

X[ H = 2 (0.158 + 0.005 (N v - 3)) ■ exp ( -™ (1 = °- M ln ^ lo) - / 2 

V 886 (l + fif(iV,-3)) 1/2 

~ 0.25 + 0.012 (^- 3) + 0.0082 ln^o (45) 

where N v is the number of massless neutrino species. This result is in a very good agreement 
with the results of the numerical calculations presented in Fig. 3 0. 

In fact, this agreement can be made even better if one notes that the formation of 4 He is 
not an instantaneous event which happens at It starts at t^' and then continues for, at least, 
50 sec (see (J3HJ)). Most of the neutrons are trapped at the end. Therefore the time delay reduces 
the amount of 4 He to Xi He ~ 0.25 exp (-50/886) ~ 0.236 that is by 1,4%. 
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figure 3 



As we see from (145)) the abundance of A He depends on the number of massless species 
N u . The presence of extra massless neutrino increases the 4 He— abundance by 1.2%. There are 
two reasons for this. Two third out of this increase is due to the dependence of the freeze-out 
concentration X* from the number of the massless species. In fact, more species one has, more 
fast universe expands at given temperature and hence the freeze-out of the neutrons occurs earlier, 
when their concentration is higher. The remaining one third has a similar nature. Namely, for 
given baryon density the nucleosynthesis happens at appropriate temperature. This temperature 
is reached earlier if there are more light species and therefore more neutrons survive until they 
will be captured. The dependence of the 4 He— abundance on the number of light species taken 
together with the results for the deuterium abundance allows us to put rather strong bounds on 
the number of unknown light particles which were relativistic at the time of nucleosynthesis. 

The helium abundance also depends on the baryon density (entropy per baryon) and 
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according to (J45|) increases by ~ 2% (numerical result ~ 2.5%) if the density is ten times higher. 
The physical origin of this dependence is very clear. In the universe with bigger concentration of 
baryons the nucleosynthesis begins earlier, at higher temperature (see (jHEI)); hence more neutrons 
survive till this time and more Helium-4 is formed. 



5 Deuterium 

To calculate the time evolution and freeze-out concentration of deuterium I will make some 
assumption which significantly simplify the consideration. The validity of these assumption can 
be checked a posteriori. 

First of all, I ignore 7 Be, 7 Li since their abundances as we will see later are always small 
compared to the abundances of 3 He and T. Second, I will assume that 3 He and T abundances 
always have quasi- equilibrium values, which are determined by condition that the "total incoming 
in appropriate reservoir flux should be equal to the outgoing flux" 11 (see Fig.l). For instance, in 
the case of 3 He it means that the amount of 3 He produced within some time interval in DD 
and Dp— reactions should be equal to the amount of 3 He destroyed during the same time in 
3 HeD and 3 Hen— reactions. This is well justified because the rate of the reactions in which 3 He 
is destroyed is high enough to take care about "quick adjustment" of 3 He— concentration to the 
change of deuterium abundance. 

The system of reservoirs with pipes, depicted in Fig.l is a " self- regulated system" with 
small adjustment time. The overall picture after the beginning of the nucleosynthesis is the 
following. When deuterium concentration reaches ~ 10~ 2 the rate of DD — reactions be- 
come comparable with the rate of the deuterium production via pn— interactions fsee (j43|0 and 
then dominates. The neutrons are taken from "np— reservoir" and send via "D— reservoir" 
along "DD and Dp— pipes" first to 113 He and T— reservoirs" and from there through n3 HeD 
and TD— pipes" to their final destination, namely, in " 4 He— reservoir" . Not all of the neutrons 
taken from "np— reservoir" reach the " 4 H e— reservoir" in the first try. Some of them "escape" 
on the way there. Namely, in "DDI and T-D-pipes" one neutron is released in the reactions 
DD — > 3 Hen, TD — > 4 Hen, comes back to "np— reservoir" and then participate in the next try 
to get " 4 ife— reservoir" . Thus after the beginning of nucleosynthesis there is a stationary flux of 
the neutrons from "pn— reservoir" to nA He— reservoir" through the system of "pipes" via interme- 
diate "D, 3 He and T— reservoirs" . The "widths of the pipes" (reaction rates) connected to 3 He 
and T— reservoirs depend on the concentration in the appropriate reservoir. For instance, the 
width of the " 3 HeD and TD— pipes is proportional, respectively, to 3 He and T— concentrations. 
If the amount of 3 He would increase/decrease compared to its quasi-equilibrium value the size 
of " 3 HeD —pipe" would be quickly adjusted (respectively increases/decreases) to bring its con- 
centration to the quasi-equilibrium value, which is determined by the condition of zero total 
flux. 

11 This condition reminds the first Kirchhoff's rule for the electric currents. 
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If the universe would not be expanding then finally nearly all neutrons would go to 
" 4 ife— reservoir" and as long as the temperature goes to zero, nothing would be left besides of 
the protons and 4 He. However, the expansion plays the role of "water-tap" for the "pipes" . At the 
moment when the reaction rates become smaller than the rate of the expansion the "water-taps" 
close and the abundances of the elements in the appropriate reservoirs freeze-out at their quasi- 
equilibrium values. The final abundances of 3 He and T are determined by deuterium freeze-out 
concentration which we have to calculate. 

Analyzing the system of kinetic equations one can find that even if 3 He and T have quasi- 
equilibrium concentrations the neutrons and deuterium concentrations not necessary satisfy the 
quasi-equilibrium conditions. Therefore, we have to derive the equations which describe the time 
dependence of the appropriate abundances X n , X D after X D reached the value ~ 1CT 2 . 

The reaction rate for the elements a, b which is equal to \ a b n a n b/ n2 B can De rewritten in 
terms of the concentrations by weight as 

' X ab X a X b (46) 



A a A b 

where A a , A b are the mass numbers of the elements a and b. The quasi-equilibrium condition for 
3 He takes then the following form: 

^DDlXp + -XopXnXp = -XiHeDXsHeXD + ~ XzH en Xzn e X n (47) 

Similar for tritium we have 

1 2 1 1 

-\ddiX d + -\3 Hen X3 He X n = -X TD X T X D (48) 
4 3 6 

I will assume that these conditions are always satisfied. 

The general kinetic equation for the rate of change of free neutrons concentration can be 
easily written if we take into account they are produced in the reactions DD — > 3 Hen and DT — > 
4 Hen and "destroyed" in the processes pn — > D'j and 3 Hen — > Tp : 

= ^DViXp + -XtdXtX d — X pn X p X n — -\3 Hen X-A He X n . (49) 

Assuming that tritium satisfies quasi-equilibrium condition (|4*K|) one can simplify this equation 
to ^ 

= -^XddX 2 d — \p„XpX n , (50) 
where as usually Add = Addi + ^dd2-) 

The appropriate equation for deuterium is derived similar by, using (}4Tj) and 
alX D 



dt 



2X pn X p X n — XddX d — 2Xrj p XrjXp. (51) 
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Expressing time through the temperature via (J19)) and substituting the numerical values 
for \ DD given by the above equations reduce to 



dX„ 



dT, 



MeV 



K(T) Vl0 {BiXn-Xl) 



(52) 



and 



dX 



= Aa- K (T) m X 2 D + R 2 X D - -R^X* (53) 

CtlMeV V ^ / 

where a = 0.86 x 10 5 and the coefficient K (T) accounts for the temperature dependence of (av) 
for DD-reactions and changes from ~ 1 to 0.5 when the temperature drops from 0.09 MeV to 
0.04 MeV. In the expressions 



R 1 = AX, 



A 



pn 



A 



(3^8) x 10" 



Ro = 2X r , 



A 



12 -- p A 



pL> 



CD 



(2.5 -=-2.3) x 10" 



(54) 



I used the experimental value for the ratio of the appropriate reaction rates; the first number 
within the brackets corresponds to the higher temperature when it changes in the interval T^eV — 
0.09^0.04. 

The system of equations (}52j) and (|53J) has an attractor solutions, which can be easily 
found if we consider Xd as a function X n (or vise versa) and rewrite the eqs. (|52 j) . (j53j) as 

dX D _ [Xp + R2X D — ^RiXn) 



dX„ 



{R x X n - X 



If X D < X n then 



X n — — (Xf) + R 2 X D ) 
Hi 



lXo 

8 X n 



+ 




(55) 



satisfies this equation up to the second order terms in Xr>/X n . The solution (|53|) is a good 
approximate solution after deuterium concentration reaches the maximal value about 10 -2 and 
begins to decrease (see Fig.2). It is valid until the moment when the neutron concentration 
drops and becomes comparable to the deuterium concentration. The solution (|B3j) describes the 
situation when the deuterium abundance satisfy the quasi-equilibrium condition. One can check 
that in this case the time derivative of the deuterium concentration in the l.h.s. of the equation 
(J53|) is small compared to every separate term in the r.h.s. of this equation. Since R2 <C R\ we 
infer from the eq. (|55|) that the deuterium and neutron concentration become comparable when 
the deuterium concentration drops to O (1) R\. Before this happens (for X D , X n > O (1) R\) the 
deuterium concentration can be expressed through the neutron concentration as 



X 



D 



R\X n 



(56) 
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Note that according to this formula the maximal possible concentration which deuterium can 
reach is Xr> ~ 10~ 2 when most the free neutrons are still not captured by light elements 
(X n ~ 0.12) . This is in complete agreement with naive estimate we got before comparing pn— 
and DD— reactions rates. When deuterium follows its quasi-equilibrium track (JoT)|) the neutrons 
concentration satisfies the equation 

dX 1 

■^~~a-K(T) Vl0 RiX n (57) 



dTueV 2 

In this case the neutrons are the "key element" which determines the quasi-equilibrium concen- 
trations of all other elements including deuterium. In other words, the neutrons regulate the 
"water-taps in the pipes connecting the reservoirs in Fig.l". The equation (|57jl starts to be ap- 
plicable at the moment when deuterium concentration grows to 10~ 2 . At this time most of the 
free neutrons are not yet trapped by the light elements and X n ~ 0.12. According to (|28p. which 
is still applicable at this time, the deuterium reaches the maximal possible concentration ~ 10~ 2 
when the temperature drops to 

TueV — 0-07 + 0.002 lnr/ 10 (58) 

After that the neutron concentration satisfies the equation (|57jl. the approximate solution of 
which is 

X n (T MeV ) ~ 0.12 exp ( l -a ■ K (T) r/ 10 i?i (T MeV - 0.07 - 0.002 In Vw )) (59) 



2 

As it follows from here, the neutron concentration decreases as the temperature drops and becomes 
equal to the deuterium concentration (~ Ri) when 

T* MeV ~ 0.07 - 0.02X^7]^ + 0.002 ln^ (60) 

It is clear that this formula is not applicable if t]iq < 0.35. In the universe with very low baryon 
density (rjio <C 1) the neutron concentration never drops below the deuterium concentration. It 
freezes-out before. In this case the nucleosynthesis is over very fast after beginning and neutron 
concentration freezes-out before the substantial part of the neutrons is converted into 4 He. After 
that the free neutrons decay. This explains why in the universe with very low baryon density 
(for instance, with t]iq ~ 10~ 2 ) the helium abundance is less than one percent (see Fig. 3). When 
I was deriving the formula (|3HJ) for 4 H e— abundance I assumed that the reactions converting the 
neutrons into A He are very efficient and able to transfer most of the available neutrons into more 
heavy elements. This means that this formula is valid only for 7710 > 0.35. From the observations 
of the luminous baryonic matter we know that 1 < rjio < 10 2 and therefore we concentrate from 
now on only on this range of parameter 7] w . All the derivations below will be done under this 
assumption. 

The neutron concentration drops to X n ~ Xd ~ R\ at the temperature T* given by ([60j). 
After that the solution (J56| is not valid anymore and the system quickly gets to another attractor 
which correspond to the quasi-equilibrium solution of the equation (JH2J)j namely, 



X n - ~]^ X D 




(61) 
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As X n /X D continues to decrease one can neglect the deviations from X n ~ X^/Ri and the 
equation (|53p for deuterium takes the following form 

d X 

l> 2a-K(T) Vl0 (X 2 D + 2R 2 X D ) (62) 



dT, 



MeV 



Now the deuterium becomes the "key element" and determines its own "fate" regulating simulta- 
neously the quasi-equilibrium concentrations of the other elements including the neutrons. Since 
R 2 practically doesn't change in the relevant temperature interval (see (j54|) ) it can be treated as 
a constant and the equation (JH2|) can be easily integrated: 



2R 2 \ (_ 2R 



X D (T)J V Xd(T 



»2 



exp 4ar] W R 2 / K (T) dT , (63) 



T 



where the temperature is expresses in MeV. When temperature goes to zero (T — > 0) the deu- 
terium concentration doesn't vanish, instead it freezes-out at X D = Ad (T — > 0) . Taking into 
account that Ad (T*) ~ R\^> R 2 and estimating the integral in ()63|) as ~ K (T*) T*, where T* 
is given by (|6U|) . we obtain the following expression for the deuterium freeze-out concentration: 

X f — (64) 

Xd ~ (exp(^o)-l) { } 

where 

A = 2aR 2 K (T*) T* (65) 

The numerical coefficient A only slightly varies with 7710- Actually when 7710 changes by two decade 
from 1 to 10 2 this coefficient increases only twice from ~ 0.1 to ~ 0.2 . The expression (|64|) fits 
very well the results of the numerical calculations presented in Fig. 3. If we want to get better 
accuracy using analytical approach we can do it taking into account the temperature dependence 
of the reaction rates. The formula ()64|) is in satisfactory agreement with the numerical results if 
we take A = 0.1 = const. At i]iq < l/A ~ 10 the good approximation for (p>4"j) is 

4^j~4x 10- Vo 1 (66) 

We see that in this range of 7710 the deuterium freeze-out concentration decreases nearly linearly 
with 7710. It is easy to understand. In this case the freeze-out concentration never drops below R 2 ~ 
10~ 5 and as it is clear from the equation (jH^j) DD— reactions always dominate over Dp— reactions 
in destroying deuterium. Therefore, the deuterium concentration freeze out when X^^X^ ~ t~ l ; 
since Add oc n B oc r/ 10 we see that in the leading order X{j should be inversely proportional to 
7710 (compare to (J34j) ). 



On the contrary, if r/ 10 > 10 the linear in Ad— term in the equation ()62j) dominates after 
Ad drops below R 2 ~ 10 -5 and after that X D cx exp (—7710 x function of T) ; hence as it follows 
from (JMj) the freeze-out concentration in this limit is 

X f D ~ 2R 2 exp {-Ar) W ) (67) 
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and decays by five order of magnitude from ~ 10~ 5 to ~ 10~ 10 when n w changes by only one 
decade from 10 to 100 (see Fig.3). In this case the freeze-out concentration is entirely determined 
by the reaction Dp — > 3 Hey which dominates over DD— deuterium destruction during the last 
stage before freeze-out. Hence in a dense universe nearly all deuterium is very efficiently burned 
down in this reaction. The deuterium abundance is very sensitive indicator of the baryon density. 
This allows us to put rather strong upper bound on r] w from observations. 



6 Helium-3 and Tritium 

Now we can calculate the freeze-out abundances of the other light elements using the quasi- 
equilibrium conditions and assuming that these conditions are still satisfied at the moment of 
deuterium freeze-out. First I consider 3 He and assume that X^ D > R 2 ~ 10~ 5 , that is, I consider 
the case of 1 < n 1Q < 10. In this case the deuterium freeze-out is determined by DD— reaction 
and happens at the time determined by condition XddX^ ~ t^ 1 . The freeze-out time for 3 He 
can be estimated requiring that the reaction 3 HeD — > 4 Hen becomes inefficient in converting the 
significant amount of 3 He into A He. This occurs when XzhcdXd ~ tl^e- Since \$HeD is in few times 
bigger than Xbd it is clear that 3 He concentration freezes-out a little bit later than the deuterium 
concentration. This means that at the moment of deuterium freeze-out the 3 HeD— reaction is 
still efficient in returning neutrons back to "np— reservoir" 12 ; hence the quasi-equilibrium solution 
for the neutrons (JoTj) derived under this assumption is still valid. Substituting X n = X^/Ri in 
(07J) we can express the quasi-equilibrium concentration of 3 He through Xp : 

^ _ 3 XddiXd + 2Xd p X p g. 

2 Xa H eD + 2 (XlHen/ Rl) X D 

After deuterium freezes-out the small leakage from "D to 3 H e— reservoir" is still able to keep 
stationary quasi-equilibrium "flow through 3 H e— reservoir". Actually in considered case the 
3 i7e-concentration is significantly smaller than X^ D and the deuterium demand needed to com- 
pensate the leakage of 3 He through " 3 HeD— pipe" is not very high. After a short time when the 
3 HeD— reaction becomes inefficient both "DDI and 3 HeD— pipes" close up simultaneously and 
the 3 He— freeze-out concentration can be obtained from (|UHj) substituting there X^ D given by (jMJ) 

Xi ~ °- 2X ° + 10 ~ 5 (69) 
3He l + 4xl0 3 x£' 1 j 

Here I first normalized all reaction rates on XsueD and then used the experimental values for the 
obtained ratios. Of course these ratios depend on the temperature. However, as one can check, 
they do not change too much , namely, not more than by factor two in the whole range of the 



12 The same is true for the tritium since Xtd — ^He 3 D- Note also that also DD— reaction at this time still 
continue actively participate in refilling the "np— reservoir" in accordance with small quasi-equilibrium neutron 
demand. 
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relevant temperatures for the two decade of baryon density. For the defmiteness I took them at 
T ~ 0.06 MeV. The obtained result is in excellent agreement with the results of the numerical 
calculations presented in Fig.3. We see that if X^ D ~ 10~ 3 then the appropriate 3 He abundance 
is ten times smaller that the deuterium abundance. It is clear from the above expressions that 
this suppression of 3 H e— abundance compared to the deuterium abundance is mostly due to the 
reactions converting 3 He into tritium, which are still very efficient at the moment of freeze-out 
because of the rather high concentration of free neutrons. In more dense universe where the 
deuterium abundance is smaller the availability of the free neutrons appropriately reduces. If, 
for instance, X B is smaller than ~ 2.5 x 10~ 4 (see the denominator in the formula (|55jl and 
([69))) then the reaction 3 Hen — > Tp does not play any significant role in determining the final 
3 He— abundance. This is why we can still use (|69|) to estimate X( H even in the universe with a 
relatively high baryon density, where X* D < 10~ 5 , although the free neutron concentration can 
significantly deviate from given by (|61|) . 

From (1691) it follows that, if XL ~ 1.2 x 10~ 5 , then the helium-3 freeze-out concentration 



^num ijuaji n luiiuws max, n yv D 
is equal to the deuterium concentration. This is in a very good agreement with the numerical 
results. In the universe with rj 10 > 10 the helium-3 is produced in the reaction Dp —> 3 Hej 
and destroyed in the process 3 HeD ^ 4 Hen. Irrespective how big is Xp these two competing 
processes give rise to the final 3 He— abundance X[ He ~ Xd p /^ H( , d ~ 10~ 5 even in the case when 
the deuterium is practically absent. Slight deviation of the numerical results from predicted here 
constant 3 He— concentration in this limit is due to weak temperature dependence of the reaction 
rates. 

Similar by, one can find that the tritium quasi-equilibrium concentration is equal to 

X T = + 2 >^£*l) Xd (to) 

Using the experimental values for the ratio of the reaction rates Xdd2/^dt — 0.01 and XsHen/Xor — 
1 we can easily understand the dependence of tritium freeze-out concentration on 7710- The formu- 
lae (jo^j) and ()70p explain why the 3 He— track in Fig. 2, in distinction from T, doesn't "repeat" the 
track of the deuterium , namely, the 3 He— concentration increases monotonically all the time, 
while T— concentration first reaches the maximum and then decreases until it gets its freeze-out 
value. 



7 Lithium- 7 and Beryllium- 7 

The quasi-equilibrium concentrations of 7 Li and 7 Be can be determined from the equations 

111 

Y^X4 HeT Xi He X T + -\v Ben X7 Be X n = -\7 Lip X7 Li X p (71) 

and 

Y^^ 4 HeHe 3 X4 He X3 He = -\l Ben Xl Be X n (72) 
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where I took into account only dominating reactions in which, respectively, 7 Li and 7 Be are 
produced and destroyed. One can check that if ry 10 > 1 then the other reactions as, for instance, 
7 LiD — > n + 2 4 ife and 7 BeD — > p + 2 4 ife for can be ignored. It immediately follows from these 
equations that 

7 Xi He f\4 HeT \ ( ^HeHe 3 v \ / VQ x 

^7 Li = — — — A T + — X 3He (73) 

LZ A p y ^7 Lip J \ A±HeT J 

The ratio of X*H^r/^ r Lip is remarkably constant within rather broad temperature interval, namely, 
it increases from ~ 2.2 x 10~ 3 to ~ 3 x 10~ 3 when the temperature drops in three times from 0.09 
MeV to 0.03 MeV. In distinction from that K% (T) = \4. He 3 He / Xi HeT varies quite significantly in 
the same temperature interval: K\ (T) ~ (5 4- 0.6) x 10~ 2 , that, is it drops nearly ten times when 
the temperature drops in three times. Substituting X± He ~ 0.25, X p ~ 0.75 in (f?3|) we obtain 

Xi Ll ~ (3 - 5) x 10" 4 (X T + K x (T) Xz He ) (74) 

To get the freeze-out concentration of 7 Li we have to substitute in this formula the appropriate 
value of Xt , K\ (T*) and X3 He at the moment when 7 Li— freeze-out. This moment can be 
evaluated analyzing the freeze-out condition for 7 Li, 7 Be— reactions. If Xd > 3 x 10~ 5 (1 < 7710 < 
5) this freeze-out occurs after the deuterium gets its final abundance. Therefore to estimate X$ L . 
in this case one can substitute the obtained above values for X[ H and in (J7I|). If 7710 = 1 
the first term there dominates. Taking into account that ~ 0.01 x X^ D ~ 4 x 10~ 6 we get 
Xi Li (r/10 = 1) — 10 -9 . As r/io increases drops and therefore Xf Li also decreases until the 
second term in (|71j) starts to dominate. After that the 7 Li abundance starts to increase. This 
increase is due to two reasons. First of all it comes from the freeze-out temperature dependence 
of K (T*) , which can be easily understood. Namely, the freeze-out temperature for 7 Li in this 
case is mostly determined by the efficiency of 7 Ben— reaction which in turns depends on the 
neutron concentration. In more dense universe the deuterium and free neutrons burned down 
more efficiently and disappear earlier (at higher temperature) than in the universe with small 
baryon density. Therefore the 7 Li concentration freezes-out at higher temperatures for which 
K\ is bigger. Second reason is the following. If 7710 > 5 the 7 Ben— reaction become inefficient 
before 3 He reaches its final freeze-out concentration. Therefore, to estimate XL . in this case we 
have to substitute in (|74*j) the actual value of Xzn e at the moment when the 7 Li— concentration 
freezes-out, which is bigger than X{ H . The numerical calculation show that after passing through 
a relatively deep minimum Xf L . (7710 = 2 4 3) ~ 10~ 10 the Lithium concentration comes back to 
~ 10~ 9 at r/10 — 10. The "trough" in X? L . — rjio dependence has a very simple explanation, 
namely, it is due to the competition of two reactions. In the universe with 7710 < 2 4 3 most of the 
Lithium-7 is produced directly as a result of A HeT— interactions. The efficiency of this process 
decreases with increase of 7710 and as 7710 becomes bigger than about 2 4-3 the reaction 7 Ben takes 
over compared to the direct 7 Li— production. In this case most of the Lithium-7 is produced via 
intermediate n7 Be— reservoir" 

The Beryllium-7 is not so important from the observational point of view. Therefore, to 
get an idea about its expected concentration we estimate the amount of 7 Be only for 1 < 77 10 < 5. 
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In this case 7 Be freezes-out after deuterium and the quasi-equilibrium solution (|61|) for the free 
neutrons is still valid at this time. Hence we get 



XL = B7%«' (T^) *4« ~0(1) 10---% (75) 
12 (X' D ) V Ar Bm / / x A 

where I used the experimental values for the ratios of the appropriate reactions; here the product 
of these ratios changes in about five times in the relevant temperature interval. At 77x0 = 1 we 
have 4~4x 10" 4 , X[ He ~ O.LX£ and correspondingly xf Be ~ 2.5 x 10~ 10 . 



8 Conclusions 



The derived above abundances of the light elements are in very good agreement with the results 
of the numerical calculations reviewed, for instance, in The analytical derivation is useful 

because it allow us to look inside "black box" (computer) and understand an interesting physics 
(for instance, attractors behavior etc.) taking place during nucleosynthesis. Without referring to 
the computer codes one can estimate the final abundances of the light elements and understand 
their dependences on the main cosmological parameters. The theoretical calculations of the 
abundances are in agreement with observations [H],[I] arid the estimates of the baryon density 
in the universe from nucleosynthesis come in impressive correspondence with CMB temperature 
fluctuation measurements. This gives strong support to the standard cosmological model. 
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